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We investigate the quantum evolution of large black holes that nucleate spontaneously in de Sitter 
space. By numerical computation in the s-wave and one-loop approximations, we verify claims that 
such black holes can initially "anti-evaporate" instead of shrink. We show, however, that this is a 
transitory effect. It is followed by an evaporating phase, which we are able to trace until the black 
holes are small enough to be treated as Schwarzschild. Under generic perturbations, the nucleated 
geometry is shown to decay into a ring of de Sitter regions connected by evaporating black holes. 
This confirms that de Sitter space is globally unstable and fragments into disconnected daughter 
universes. 



I. INTRODUCTION AND SUMMARY 
A. Anti-evaporation and proliferation 

Schwarzschild-de Sitter black holes have unusual quantum properties and instabilities [Q-D . They are of cosmological 
interest because they can be produced during inflation |9j . They are also theoretically significant because they change 
the global structure of de Sitter space fundamentally. 

Unlike their asymptotically flat cousins, Schwarzschild-de Sitter black holes are surrounded by a cosmological 
horizon. This limits their size and puts them in a thermal bath. Their temperature is always larger than that of the 
cosmological horizon Nevertheless, Hawking and one of the authors argued in Ref. Q that some black holes 

accrete so much quantum radiation that they will grow, or 'anti-evaporate,' instead of evaporating (see also 

In the maximal solution, the black hole and cosmological horizon are of equal size, and the spatial geometry will 
be 5^ X S*^, with constant two-sphere radius. This geometry can nucleate semiclassically in de Sitter space through a 
gravitational tunneling process. Its evolution is unstable to the formation of n de Sitter regions, distributed around the 
5'"'^ and connected by n black hole throats 0. This configuration may be visualized as a 'doughnut' with n 'wobbles,' 
whose minima and maxima represent the black hole and cosmological horizons, respectively. For n = 1, this reduces 
to the ordinary, submaximal Schwarzschild-de Sitter solution. For higher n, it resembles a ring-like sequence of n 
Schwarzschild-de Sitter solutions. As the de Sitter regions form between the cosmological horizons, the black holes 
are expected to evaporate and shrink. When they finally disappear, the doughnut is pinched in n places, leaving n 
disconnected pieces. This corresponds to the fragmentation of space into n large daughter universes [Q. 

Thus de Sitter space is globally unstable. Locally, however, the daughter universes are indistinguishable from 
de Sitter. Therefore, they will harbour more maximal black holes, whose evaporation can lead to further fragmentation. 
Iteratively, an unbounded number of disconnected de Sitter universes is produced: de Sitter space proliferates. This 
drastic conclusion depends crucially on the nonlinear evolution of the instabilities of the maximal Schwarzschild- 
de Sitter solution (also referred to as the 'Nariai' solution below). 



B. Analytical shortcomings 

In Refs. and 0,^, the anti-evaporation and proliferation effects were derived analytically, using a number of 
approximations. The stability of the Nariai geometry was studied using linear perturbation theory. The evolution of 
the horizons was found only for very early times, as a power series in the time variable. This led to the discovery of 
anti-evaporation for a class of perturbations. But a power series can prove only initial anti-evaporation. It cannot 
answer an important question; Will such black holes continue to grow, asymptotically approaching the maximal 
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(Nariai) solution? Or will their growth eventually stall, and evaporation set in after all? This problem is highlighted 
by recent claims that anti-evaporation may permanently stabilize black holes nucleated during inflatiorj^ [01 . 

At least for Euclidean boundary conditions, it was shown that the black holes do evaporate at late times H. In fact, 
arguments might be made that the evaporating mode is attractive at late times for all initial conditions In any 
case, however, the intermediate period containing the putative transition between anti-evaporation and evaporation 
is beyond the reach of the analytical methods employed. Can there be a smooth transition, and how does it proceed? 
To understand this crucial phase of the black hole evolution, it should be explored numerically. 

If anti-evaporation is only temporary, then how can we be sure that the evaporating phase, once entered, continues 
until the black hole disappears? This question is particularly important in the context of multiple black hole formation 
and proliferation The fragmentation of space is possible only if all black holes become small and disappear like 

Schwarzschild black holes. 

To show that black holes really evaporate completely, it is necessary to trace the backreaction at least until they 
are significantly smaller (perhaps by a factor of 10) than the cosmological horizon. Then they can be assumed to 
behave like Schwarzschild black holes, which presumably radiate all of their mass away. But linear perturbation theory 
requires that the small two-spheres inside the black hole and the large two-spheres between the cosmological horizons 
still be of nearly equal size. Thus, the analytic 'late-time' solution, showing black hole evaporation, is not really valid 
for arbitrarily late times. A full, non-linear numerical calculation is needed to interpolate to the Schwarzschild regime. 

C. Numerical analysis 



In Sec. [II we consider n — 1 perturbations of the Nariai geometry, which can be parametrized by their initial 
amplitude, e, and phase, Depending on the black holes undergo initial periods of anti-evaporation or evaporation. 
We demonstrate the transitory nature of anti-evaporation for neutral black holes by showing that evaporation sets 
in at sufficiently late times, independently of ??. The evaporation rate is compared to an approximate solution of the 
linear equations in order to verify that this solution is an attractor. Beyond the linear regime, the black holes are 
found to shrink to a size much smaller than the cosmological horizon before our computation breaks down. Then they 
live in approximately flat space, and should simply continue to evaporate until they disappear. 



In Sec. IV we show that maximal black holes become stable if their charge exceeds a critical value, in agreement 
with thermodynamic expectations. This confirms analytic arguments given in Ref. In particular, we show that 
supercritical black holes, independently of their initial behavior, always end up anti-evaporating. This is the condition 
for a novel type of daughter-universe production proposed in Ref. 

In Sec. we show that a higher mode (n > 1) perturbation behaves as predicted by analytical arguments: First 
it oscillates, then it freezes out, forming black hole interiors and growing de Sitter regions. Depending on initial 
conditions, a period of anti-evaporation may follow. Finally, the black holes evaporate and become much smaller than 
the cosmological horizons. 

II. ACTION AND PERTURBATIONS 
A. Including bacli-reaction 

The four-dimensional Lorentzian Einstein- Hilbert action with a cosmological constant. A, and a Maxwell field, F^i,, 
is given by: 

S=^jd'x (-gi^)i/^ [i?^^ - 2A - F,,F^^] , (1) 

where BF^ and g^^ are the four-dimensional Ricci scalar and metric determinant. 

Restricting to spherically symmetric fields and quantum fluctuations, the metric may be written as 

ds^ = e^P {-df + dx^) + e-^*dn^, (2) 



^We point out, however, that even stable black holes nucleated semiclassically during inflation would_be diluted by the 
cosmological expansion. Typically, they wiU lie outside the current horizon and will not be observable 
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where x is the coordinate on the , with period 27r. Using this ansatz, and the on-shell condition for magnetic fields, 



F^.F^^" = 2Q2e4*, (3) 
the angular coordinates and the Maxwell field can be integrated out in Eq. (|^), which reduces the action to 

S^-^ I (fx (-5)i/2e-2^ \R + 2(V0)2 + 2e2<^ - 2A - 2Q^e^n , (4) 

lOTT J 

In order to include back-reaction effects, a Polyakov term, which arises in the one-loop effective action of a two- 
dimensional scalar field, will be included. (It would be preferable to work with dilaton-coupled scalars [pl]-p2[, or 
even better, with a four-dimensional effective action reduced to two-dimensions. This would complicate the numerical 
computation enormously and will not be attempted here. For small perturbations, it has been shown that extra 
terms from dilaton-coupling do not affect results 1^,^ . Thus one would not expect qualitative changes even for black 
holes noticeably smaller than the cosmological horizon. In order to corroborate our results for nearly-Schwarzschild 
black holes, however, a full four-dimensional treatment would be desirable in future work.) In the large N limit, 
the contribution from the quantum fluctuations of the scalars dominates over that from the metric fluctuations. In 
order for quantum corrections to be small, one should take iVA ^ 1. One can obtain a local form of this action 
by introducing an independent scalar field Z which mimics the trace anomaly [||. The on-shell equivalence of the 
equations of motion can be seen by choosing a conformal gauge for the two-dimensional metric, as in Eq. (2). Thus 
the action of the one-loop model will be given by: 



-2(j> 



N 



-Z 



i?--(VZf 



2 + 2e-^'^ (V(/))^ - 2e-2'^A - 2Q 



(5) 



B. Equations of motion 

Differentiation with respect to t {x) will be denoted by an overdot (a prime). For any functions / and g, define: 



a/a<? = -/g + /v, = (6) 

SfSg = fg + f'g', 6'g = g + g". (7) 
Variation with respect to p, (f> and Z yields the following equations of motion: 

- + 2{dcj)f + -e^^d^Z + e^P+^^ (Ae-^^ + Q^e^^ - l) = 0; (8) 
6 

d^p - d^(f> + {d(t>f + Ke^P - Q2e2p+4^ = 0; (9) 

d^Z - 2d^p = 0. (10) 

The constraint equations are: 

((^20 - 25(1) Sp) - (6^)'^ = ^e^'^ [{SZf + 2S'^Z - AdZSp] ; (II) 

- p(t)' - p'4^ - W = ^e^^ ZZ' + 2Z' -2 i^pZ' -f p'z) . (12) 



From Eq. ([lO|), it follows that Z = 2p + rj, where rj satisfies d'^g = 0. The remaining freedom in g can be used to 
satisfy the constraint equations for any choice of p, p, (/) and ((> on an initial spacelike section B. 
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C. Metric and horizon perturbation 



Maximal (Nariai) black holes nucleate semiclassically in de Sitter space. This process is mediated by gravitational 
instantons and has been described in detail in Refs. Here we are interested not in the nucleation, but in the 

further evolution of the Nariai solution. Its metric is given by 



.2p _ 



1 1 

A cos^< ' 



where B is given by the cubic equation 



B 



1-Q-'B[1 + ^B 



N 

1 B 
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A; 



the physical solution is the one that limits to i? = A for N = Q ~ 0). A is given by 

A = A-Q^B^. 



(13) 



(14) 



(15) 



Quantum fluctuations will perturb this solution, so that the two-sphere radius, e will vary slightly along the 
one-sphere coordinate, x. Decomposition into Fourier modes on the yields the perturbation ansatz 



p20 _ 



A. 



l + 2e^(a„(i) cos nx + dn (t) sin nx) 



(16) 



where e is taken to be small. This will be referred to as the metric perturbation, characterized by the (7„, (t„ at the 
time t — 0. 

More generally, one should also consider perturbations of the time derivative of the two-sphere radius, expressed 
by &n ■ As in Ref . B , we parameterize the initial conditions for the perturbation as 



(T„(0) = sinz? 



,(0) = cos-d 



(17) 



where i9 represents the phase of the initial perturbation. 

Based on linear perturbation theory, one would expect such perturbations to lead to a classical, as well as a quantum 
instabihty [|J|J|,|,|,|. Classically, the regions on the where the two-spheres are smaller than the Nariai value 
should collapse to form black hole interiors. The larger two-spheres, on the other hand, should grow exponentially, 
developing into asymptotic de Sitter regions. If the first mode dominates, this simply leads to a nearly-maximal 
Reissner-Nordstrom-de Sitter solution. But if higher modes are strongly excited, a whole sequence of Reissner- 
Nordstrom-de Sitter solutions can develop around the same one-sphere. This may be thought of as a necklace of 
de Sitter regions, strung together by black hole throats. 

The expected quantum instability is the evaporation of these black holes. Actually, we will show below that black 
holes of sufficient charge are stable. But all other black holes would be expected to evaporate and get smaller, since 
their temperature is higher than that of the surrounding cosmological horizon. 

The black hole evolution can be calculated numerically using Eqs. (||)~(|l^). It is important to stress that the black 
hole and cosmological horizons do not in general correspond to the minimal and maximal two-spheres along the S^, 
although such a slicing can always be found. In general, one must first find the positions of the horizons on the 
one-sphere; this can be done by finding the points where the gradient of the two-sphere size is null Q. The two-sphere 
sizes at those locations give the size of the black hole. The black hole evolution can be monitored by following the 
horizon location and plotting the horizon size vs. time. It will be convenient to define the horizon perturbation 6, for 
a black hole located at Xh, by 



Mt)-^ = e2'^[*^^'^(*)l =B[1 + 2eS{t)] . 



(18) 



Thus S corresponds to the fractional difference between the current black hole size and the size of a maximal black 
hole of equal charge. Evaporation corresponds to increasing values of 5. 

Because the expands exponentially, the black hole and cosmological horizons of a Reissner-Nordstrom-de Sitter 
solution appear to move ever more closely together in comoving coordinates (see Fig. |^). 
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FIG. 1. Location of the cosmological and black hole horizons on the coordinate for the perturbation mode n — 1. 



D. Numerical technique 



The equations of motion - (|10|) were solved numerically using the standard method of characteristics for second 
order quasi-linear hyperbohc equations ]2^ . For this purpose, the auxihary field Z can be eliminated by inserting 



Eq. ( |10| ) into Eq. (|9|) and the remaining equations can be rewritten in manifestly hyperbolic form: 



(19) 



(20) 



where 



C 



(21) 



The charateristic curves, along which the equations reduce to ordinary differential equations, are x± — x :kt. Given 
the solution at two neighboring points on non- identical characteristics and approximating dx± by Ax±, the solution 
at the intersection of the characteristic curves going through these points can be found by iteration. Assigning the 
initial conditions Eqs. ( |l3| ) - ( |l6| ) on a t = hypersurface thus allows the advancement of the solution along the 
a;±-grid, employing periodic boundary conditions at a; = and x = 2tt. 

For most of our computations, we used a resolution of 10000 grid points on x-hypersurfaces. The accuracy of the 
results was verified by comparison with perturbative solutions and with results obtained from an independent 
second-order finite difference code for the same equations. With the exception of the very last time step, the closest 
to the conformal time coordinate singularity, excellent agreement of all solutions was found. 



III. ANTI-EVAPORATION AND TURNAROUND 



Anti-evaporation was first found in Ref. Q for the § — -k /2, n — 1 perturbation of the maximal (Nariai) neutral black 
hole solution in de Sitter space. A power series approximation showed that black holes formed by this perturbation 
will grow initially. A numerical calculation, however, allows us to probe the black hole evolution beyond the range of 
validity of the power series. Interestingly, this proves the anti-evaporation effect to be transitory. This is demonstrated 
by the long-dashed line in Fig. ^j. The horizon perturbation first decreases quadratically, signaling anti-evaporation, 
in quantitative agreement with the result of Ref. B. At a time t k, .2, however, it turns around and starts to 
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increase. This means that the black hole eventually stops growing. Instead, it starts to shrink in size, corresponding 
to evaporation at late times. 

For other initial phases, the behavior of the black hole can be quite complicated, going through various periods of 
evaporation and anti-evaporation (Fig. H). The important result is, however, that it always ends up evaporating at 
late times, as conjectured in Ref. |8|. The numerical result for the final evaporation rate agrees with the asymptotic 
solutions given there to within better than one percent: 



6 cx (cosi)i~'=+. 



where C-|_ is the larger root of 



c(c+l) = 2\/l + ^7VA+iA^2A2. 



(22) 



(23) 



An initial phase just under = 7r/2 corresponds to a large time-derivative of the metric perturbation directed 
opposite to the perturbation. In this case, the metric overshoots the Nariai value and a cosmological horizon forms. 
The corresponding black hole is located where the cosmological horizon would usually sit. Its evaporation is mirrored 
by the growth of the cosmological horizon. In Fig. ^ this cosmological horizon shows up for = 57r/6 as a negative 
value of S which increases in magnitude. 
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FIG. 2. Evolution of the black hole horizon pertur- 
bation 5 for different initial values of ■£>: — (solid), 
i9 = 7r/6 (dotted), = 7r/3 (short dashed), ■& = -k 12 
(daBhed-single dotted), i9 = 2n/3 (dashed-triple dotted), 
■d = 57r/6 (long dashed) for A = 0.1, = 5, and w = 0. 
The behavior is independent of the initial amplitude e as 
long as e ^ 1; all plots are for e = .0005 unless otherwise 
indicated. 
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FIG. 3. Numerically one can show that perturbed 
Nariai black holes evaporate until are significantly smaller 
than the cosmological horizon. Then they can simply be 
viewed as living in asymptotically flat space, so evapo- 
ration will continue. Model parameters here are e — 0.5, 
1? = 7r/3. For this run, an RST-type counterterm [24] was 
included in the effective action in order to ensure that the 
black hole singularity resides at r = 0. 
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How can we be sure that the evaporation seen in the later stages of Fig. g is permanent? By the nature of the 
conformal time i, an infinite amount of proper time is contained in the final step of the numerical calculation. No 
matter how refined the step size, this region cannot be resolved, ft is sufficient, however, to trace the evolution until 
the black hole size differs significantly from the size of the cosmological horizon. This is precisely the condition for 
non-linearity, so it cannot be verified analytically using the power series for S. Our numerical calculation, however, 
can reach this regime (see Fig. Once it is so small, the black hole will be much hotter than the cosmological 
radiation, and will no longer be influenced by it. It will behave like a Schwarzschild black hole, so we can trust that 
it will continue to evaporate. 
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FIG. 4. Evolution of the horizon perturbation 5 
for black holes with charges Q^A — 0.195 (solid), 
Q^A = 0.215 (dotted), Q^A = 0.235 (dashed), 
Q^A = 0.245 (dashed-single dotted), Q'^A = 0.255 
(dashed-triple dotted), and Q^A — 0.265 (long dashed). 
Model parameters were A = 0.1, N = 5, and ^ — 0. 
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FIG. 5. Evolution of the horizon perturbation S for 
slightly supercritically charged (Q^A = 0.265) black holes 
with different initial values of 1?: 1? = (solid), ■& = n/6 
(dotted), ■& — 7r/3 (dashed), i9 = 7r/2 (dashed-single dot- 
ted), i9 = 27r/3 (dashed-triple dotted), i9 = 57r/6 (long 
dashed); model parameter values as in Fig. ^. In this 
case the Nariai geometry is an attractor. 



IV. STABILITY OF CHARGED BLACK HOLES 



Nariai solutions can be electrically or magnetically charged. In this case, their spacelike sections will still be given 
by a direct product of x S*^ (with constant radius), but the 5*^ and S"^ radii will not be equal. The flux 
runs around the . The metric is classically unstable to small perturbations. Reissner-Nordstrom-de Sitter black 
holes form where the two-sphere size is smaller, and asymptotic de Sitter regions develop where the two-spheres are 
larger |,|. 

These black holes are of nearly maximal size. For a small charge, one would expect their initial behavior to be 
similar to the neutral Schwarzschild-dc Sitter black holes. When they are highly charged, however, they become 
nearly extremal. The temperature of a black hole decreases as it approaches extremality. Thermodynamic arguments 
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were given in Ref. which showed that there is a critical value of the charge Qc = 3/(16A) = 3/4(5^^^^^^. Roughly 
speaking, when a subcritical Nariai solution is perturbed, the temperature of the black hole will be larger than that 
of the cosmological horizon because it is smaller. When a supercritial Nariai solution is perturbed, the black hole will 
be colder than the cosmological horizon, even though it is smaller, simply because its mass is already very close to 
the mass of the extremal, zero-temperature solution. 

Thus one would expect subcritically charged black holes to evaporate at late times. Supercritically charged black 
holes, on the other hand, should absorb quantum radiation from the cosmological horizon and grow, their size 
approaching the Nariai value asymptotically. This is confirmed numerically, as shown in Figs. ^ and |5[ The value 
of the critical charge is confirmed quantitatively. This is a non-trivial check that this simple, two-dimensional model 
refiects the thermodynamic properties of Reissner-Nordstrom-de Sitter black holes accurately. 

The anti-evaporation of supercritially charged black holes has an important consequence for the quantum global 
structure of de Sitter space. It gives rise to a new type of proliferation effect, as pointed out in Ref. ||^. Because 
the two-sphere size is forever nearly constant in the Hubble-size region between the two horizons, a small quantum 
fluctuation can easily invert the role of the horizons. In other words, it can increase the two-sphere size on the black 
hole horizon, turning it into a cosmological horizon, and vice- versa. This amounts to the insertion of a new black hole 
'bead' into the 'necklace.' Such processes repeat endlessly, so that an unbounded number of causally disconnected 
de Sitter regions develop on the same S^. 



V. HIGHER MODES AND PROLIFERATION 



We now return to uncharged solutions, but consider modes with n > 1. Such perturbations lead to a spatial 
geometry that can be described as a 'doughnut with n wobbles.' Like for n — I, the metric perturbation is classically 
unstable. The mode will oscillate until the expansion has stretched it enough to leave the horizon. Then it 
will freeze out, and grow exponentially. This was shown in a linear approximation in Ref. Q], and is demonstrated 
numerically in Fig. |6| for n — 3. 




FIG. 6. Metric perturbation evolution for n = 3 as function of x and t. Model parameter values as in Fig. |2[ 

The evolution of the horizon perturbation for n = 3 is shown in Fig. |^. While the metric perturbation oscillates, the 
(V0)^ — surfaces move rapidly and 'cross over' (see Fig. ||). They only represent black hole horizons after the metric 
perturbation freezes out. Fig. |^ shows that in the n = 3, = 7r/2 case considered there, the metric perturbation 
freezes out while it of the opposite sign compared to its initial value. This means that cosmological horizons develop 
from the initial minimal two-spheres. This is reflected in the figure in the negative values for the 'black hole horizon' 
perturbation. Note that the absolute value is increasing at late times, signalling evaporation. 
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FIG. 7. Evolution of the black hole horizon perturba- 
tion 5 for the perturbation mode n = 3. Model parameter 
values as in Fig. |2| 
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FIG. 8. Location of the cosmological and black hole 
horizons on the coordinate for the perturbation mode 
n = 3. 



ACKNOWLEDGMENTS 

We thank the Aspen Center for Physics where this investigation was initiated, and the German American Academic 
Council (GAAC) for financial support. 



[1] 

[2] 

[3] 
[4] 

[5] 
[6] 

[7] 
[8] 
[9] 
[10] 

[11] 
[12] 

[13] 
[14] 
[15] 
[16] 



G. W. Gibbons and S. W. Hawking: Cosmological event horizons, thermodynamics, and particle creation. Phys. Rev. D 
15, 2738 (1977). 

P. Ginsparg and M. J. Perry; Semiclassical perdurance of de Sitter space. Nucl. Phys. B222, 245 (1983). 
R. Bousso and S. W. Hawking: Pair creation of black holes during inflation. Phys. Rev. D 54, 6312 (1996), gr-qc/9606052 . 
R. Bousso and S. W. Hawking: (Anti-) evaporation of Schwarzschild-de Sitter black holes. Phys. Rev. D 57, 2436 (1998), 



hep-th/9709224 



S. Nojiri and S. D. Odintsov: Quantum evolution of Schwarzschild-de Sitter (Nariai) black holes, |hep-th/9804033 
E. Elizalde, S. Nojiri and S. D. Odintsov: Possible quantum instability of pri mordial black ho les, hop-th/9901026 
R 
R 
R 
R 
E. 



hep-th/9902183 



Bousso: Proliferation of de Sitter space. Phys. Rev. D 58, 083511 (1998), |hep-th/9805081 
Bousso: Quantum global structure of de Sitter space, accepted for publication in Phys. Rev. D 

Bousso and S. W. Hawking: The probability for primordial black holes. Phys. Rev. D 52, 5659 (1 995), gr-qc/9506047 
Bousso: Pair creation of dilaton black holes in extended inflation. Phys. Rev. D 55, 4889 (1997), gr-qc/9610040. 
Elizalde 



S. Naftulin and S. D. Odintsov: Covariant effective action and one loop renormalization of 2-D dilaton gravity 
with fermionic matter. Phys. Rev. D 49, 2852 (1994), |hep-th/930802C| . 

V. Mukhanov, A . Wipf and A. Zelnikov: On 4-D Hawking radiation from effective action. Phys. Lett. B332, 283 (1994), 
hep-th/9403"0l| . 

T. Chiba and M. Siino: Disappearance of black hole criticality in semiclassical general relativity. Mod. Phys. Lett. A12, 
709 (1997). 
R. Bousso and S 
hep-th/970523^ . 



W. Hawking: Trace anomaly of dilaton- coupled scalars in two dimensions. Phys. Rev. D 56, 7788 (1997), 



S. Nojiri and S. D. Odintsov: Trace anomaly and nonlocal effective action for 2-d conformally invariant scalar interacting 
with dilaton. Mod. Phys. Lett. A12, 2083 (1997), |hep-th/9706009| . 

S. Nojiri and S. D. Odintsov: Trace anomaly induced effective action for 2-D and 4-D dilaton coupled scalars. Phys. Rev. 
D 57, 2363 (1998), 



hep-th/9706143 



9 



[17] 
[18] 

[19] 
[20] 
[21] 

[22] 
[23] 
[24] 



S. Ichinose: Weyl anomaly of 2-D dilaton- scalar gravity and hermiticity of system operator, hep- th/9707024 
W. Kummer, H. Liebl and D. V. Vassilevich: Hawkin g radiation for n onminimally coupled matter from generalized 2-D 
black hole models. Mod. Phys. Lett. A12, 2683 (1997), |hep-th/970704l[ 

J. S. Dowker: Conformal anomaly in 2-d dilaton scalar theory. Class. Quant. Grav. 15, 1881 (1998), [hep-th/9802029 . 



;),Fe 

2, ^ 



Y. Gusev and A. Zelnikov: Two-dimensional effective action for matter fields coupled to the dilaton, pep-th 9910198 . 

F. C. Lombardo, F.D. Mazzitelli and J.G. Russo: Energy-momentum tensor for scalar fields coupled to the dilaton m two 
dimensions. Phys. Rev D59, 064007 (1999). 

R. Balbinot and A. Fabbri: Hawking radiation by effective two-dimensional theories. Phys. Rev. D59, 044031 (1999). 

G. D. Smith: Numerical Solution of Partial Differential Equations, Oxford: Clarendon Press (1985). 

J. G. Russo, L. Susskind and L. Thorlacius: Black hole evaporation in 1 + 1 dimensions. Phys. Lett. B292, 13 (1992), 



hep-th/9201074 



10 



